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ABSTRACT 

We explore the cosmic evolution of massive black hole (MBH) seeds forming within 
'quasistars' (QSs), accreting black holes embedded within massive hydrostatic gaseous 
envelopes. These structures could form if the infall of gas into the center of a halo 
exceeds about 1 M© yr -1 . The collapsing gas traps its own radiation and forms a 
radiation pressure-supported supermassive star. When the core of the supermassive 
star collapses, the resulting system becomes a quasistar. We use a merger-tree approach 
to estimate the rate at which supermassive stars might form as a function of redshift, 
and the statistical properties of the resulting QS and seed black hole populations. We 
relate the triggering of runaway infall to major mergers of gas-rich galaxies, and to a 
threshold for global gravitational instability, which we link to the angular momentum 
of the host. This is the main parameter of our models. Once infall is triggered, its 
rate is determined by the halo potential; the properties of the resulting supermassive 
star, QS and seed black hole depend on this rate. After the epoch of QSs, we model 
the growth of MBHs within their hosts in a merger-driven accretion scenario. We 
compare MBH seeds grown inside quasistars to a seed model that derives from the 
remnants of the hrst metal-free stars, and also study the case in which both channels of 
MBH formation operate simultaneously. We find that a limited range of supermassive 
star/QS/MBH formation efficiencies exists that allows one to reproduce observational 
constraints. Our models match the density of z — 6 quasars, the cumulative mass 
density accreted onto MBHs (according to Soltan's argument), and the current mass 
density of MBHs. The mass function of QSs peaks at Mqs — 10 6 M©, and we calculate 
the number counts for the JWST 2 — 10 ftm band. We find that JWST could detect 
up to several QSs per field at z ~ 5 — 10. 

Key words: black hole physics — infrared: stars — cosmology: theory — galaxies: 
formation — galaxies: nuclei — quasars: general 



1 INTRODUCTION 

While there is ample evidence that supermassive black holes 
populate the nuclei of most large galaxies and that some 
black holes w ith masses exceeding 10 9 Mq formed as early 
as z > 6 fe.g.,rFarJl200ll:lB"arth et alJl2003:lDiorgovski et all 



|200S| ; IWillott et all [20091 : 1 Jiang et ail 120091 ), there is little 
consensus as to the progenit ors of these h oles. Two schools of 
thought have persisted since lReesl (1978) first devised a flow 
chart outlining possible routes of massive black hole (MBEQ) 
formation. According to one line of argument, supermassive 
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1 We refer here generically to MBHs when the hole mass is above 
the limit for black hole formation in today's stars, ~ 50 Mq. This 
definition comprises both seed black holes and supermassive black 
holes in galaxies and quasars. 



black holes grew from the remnants of an early population of 
massive stars, the so-called Population III (Pop III), which 
is believed to have formed in pregalactic minihalos at z > 
20. According to the other, the precursors of supermassive 
black holes could have formed by the 'direct collapse' of large 
amounts of gas in much larger halos at later times. 

Each scenario has both positive and negative at- 
tributes. Although Pop III remnants were unlikely to have 
been more massive than a few hundred Mq each, they 
would have formed relatively early and thus their growth 
process would have had a considerable head start. They 
could have congregated and merged in the cores of merging 
minihalos, while simultaneously growing by accretion. 
Models for the gro wth of superma s sive bl ack holes from 
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supermassive black holes, but have rather more difficulty 
in producing enough 10 9 Mp holes at z > 6 to explain the 
earliest known quasars (IVolonteri fc Reed 120051 : IShapirol 
120051 : IVolonteri fc Reesll2006r i. Additional worries about the 
Pop III scenario include the possibility that too many of the 
remnants would have been ejected from the cores of merging 
halos and that their accretion rates would be depressed 
by the shallow potential wells of the host mini-halos and 
heating of the ambient g as by stellar radiation and winds 
|Milosavlievic et al.ll2009l . and references therein). 

Seed black holes that formed by direct collapse 
(e.g., at z < 15) would have had less time to grow, 
but this would have been partly compensated by their 
larger initial masses. Various direct co llapse models 
for s ee d formation have been pr oposed (lLoeb fc R asio 



19941; lEisenstein fc Loebl 1 19951; lHaehnelt et al] 1998 
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2004 



2003; Koushiappas ct al 
Begelman et al.ll2006l ; lLodato fc Nataraianll2006l ; IBegelmanl 
2010h , but so far only limited attempts have been made 



to place these models into the context of st ructure for- 
mation theories (e.g., IVolonteri et alj 120081 ; lLagos et al.l 
120081 ; IVan Wassenhove et al.ll2010h . Because direct collapse 
models draw only indirectly on star formation lore, there is 
much less consensus about the initial conditions for direct 
collapse and the details of how it might have occurred. 

In a series of recent papers, we have described a se- 
quence of events that we believe repre sents a plausible route 
to MBH formation via direct collapse fecgclma n et al . 2006, 
120081 ; [Begelman! l20ld) . In this picture, the main triggering 
event is the infall of gas into the center of a halo at a rate 
exceeding about 1 Mq yr _1 . Such large rates of infall are 
possible in halos with virial temperatures in excess of about 
10 4 K, which only become common at z < 10—15. Questions 
remain about the ability of gas to accumulate at such a high 
rate without too much of it fragmenting into stars, but as 
we argue below, recent simulations as well as analytic cal- 
culations suggest that the importance of fragmentation may 
have been severely overestimated in past work. The collaps- 
ing gas traps its own radiation and forms a quasistatic, ra- 
diation pressure-supported supermassive star, which burns 
hydrogen for about a million years while growing to a mass 
> 10 6 Mo (Begelman 2010; note that the earlier claim in 
Begelman et al. 2006 that H-burning is unable to postpone 
collapse is erroneous and is corrected in the later paper). 

Because of rotation, the black hole that forms initially 
probably comprises only a small fraction of the core, with 
a mass < 1O 3 M0, but it grows rapidly at a rate set by 
the Eddington limit for the massive gaseous envelope. A 
novel feature of our model is the prediction that the envelope 
swells by a factor of > 100 in radius as it absorbs the energy 
liberated by black hole growth. The resulting object, which 
we have dubbed a 'quasistar', (QS) resembles a red giant 
with a luminosity comparable to a Seyfert nucleus. As the 
black hole grows inside it, its photosphere expands and cools 
until it hits a minimum temperature associated with the 
Hayashi track, at which point it disperses, leaving behind 
the naked seed black hole. 

In this paper, we use a merger-tree approach to esti- 
mate the rate at which supermassive stars might form as 
a function of redshift, and the statistical properties of the 
resulting QSs and seed black holes. In Section 2 we sum- 
marize the properties of supermassive stars, QSs and seed 



black holes, as a function of time and gas infall rate. The ex- 
istence of runaway infall depends on a threshold for global 
gravitational instability, which we model as a function of 
the gas mass and angular momentum following each halo 
merger — we discuss these criteria in Section 3. Once infall 
is triggered, its rate is determined by the halo potential; the 
properties of the resulting supermassive star, QS and seed 
black hole depend on this rate. 

To relate the properties of the seed black holes to the 
observable distributions of active and quiescent black holes 
at different redshifts, we apply a simple set of rules for their 
subsequent growth by accretion and mergers. These are de- 
scribed in Section 4, and in Section 5 we present our results. 
We compare our direct collapse seed model to a pure Pop 
III seed model, and also study the case in which both chan- 
nels of MBH formation operate simultaneously. We summa- 
rize our conclusions and discuss the prospects of testing this 
model observationally in Section 6. 



2 FORMATION OF SUPERMASSIVE STARS, 
QUASISTARS AND SEED BLACK HOLES 

Theoretical aspects of the formation and evolution of super- 
massive stars, QSs, and the resulting seed bl ack holes are dis- 
cussed in Begelman e t al. 20 06 fBVR2006'). lBegelman et al.1 
|2008l ) and IBegelmanl (|2010h . Here we summarize the formu- 
lae that we use as inputs for merger-tree calculations pre- 
sented below, with emphasis on the results from Begelman 
et al. (2008) and Begelman (2010). 



2.1 Supermassive stars 

A supermassive star in the center of a forming galaxy is 
unlikely to evolve at constant mass, but rather accumu- 
lates mass due to infall. We adopt a constant accretion rate 
M = m Mq yr _1 . For infall rates larger than about 1 Mq 
yr _1 , the star never relaxes thermally to a fully convective 
state during its "main sequence" (hydrogen-burning) life- 
time; instead it develops a convective core (with polytropic 
index n = 3) surrounded by a convectively stable envelope 
(Begelman 2010). The main sequence lifetime is then 



t MS ~ lAm- 1/2 T-l Myr, 



(1) 



where 10 8 T Ci 8 K is the central core temperature, which de- 
pends weakly on the mass and CNO abundance as shown in 
Fig. 5 of Begelman (2010). For the calculations below, we 
adopt a fiducial value of T Cj 8 = 1.4, which corresponds to 
a CNO abundance by mass Z C no ~ 10" 6 (Z/Z Q ~ 10" 4 ). 
Then, £ms ~ m _1//2 Myr and the star's mass at the end of 
hydrogen burning is 



10 6 m 1/2 Mr. 



(2) 



A partially convective supermassive star has a radius 
R, ~ 3.4 x 10 13 m cm (for electron scattering opacity 
k = 0.34 cm 2 g _1 ), independent of mass, and an effective 
temperature 



T. = 1.2 x 10 5 m _3/8 K 



(3) 



at the end of hydrogen burning. Although the values of 
T* and i?» suggest that supermassive stars may be strong 
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sources of far ultraviolet radiation, we cannot confidently 
calculate the ionizing flux from a rapidly growing supermas- 
sive star because the infalling matter is very optically thick. 
The opacity is strongly scattering-dominated so there will 
be a color correction to the effective temperature that hard- 
ens the spectrum, but if the infalling matter obscures the 
supermassive star then the true photosphere could be well 
outside the star's nominal radius, leading to softer radiation. 
By using a blackbody with temperature T* and radius R, 
to estimate ionizing fluxes, we are effectively assuming that 
the infalling matter is distributed highly anisotropically in 
solid angle (e.g., due to rotation). 

Because they are strongly dominated by radiation pres- 
sure, nonrotating supermassive stars have very small bind- 
ing energies and can b e destabilized by tiny corrections d ue 
to general r elativity (|Hovle fc Fowlerlll963l ; llbenl Il963al lbl; 
lFowlerlll96a ). However, a relatively small amount of rotation 
(compared to breakup) can increase the binding energy dra- 
matically, stabilizing stars as massive as ~ 10 8 Mq (Fowler 
1966). The gas falling onto supermassive stars is expected 
to have substantial, though sub-Keplerian, rotation, and we 
will conservatively assume that supermassive stars can ex- 
ist up to masses of 10 7 Mq , above which they undergo core 
collapse even if hydrogen is not exhausted. 

2.2 Quasistars and seed black holes 

It is probable that only a tiny fraction of the mass in the core 
of a supermassive star has low enough angular momentum 
to collapse directly to a black hole. Any additional growth of 
the hole must be accompanied by the outward transport of 
angular momentum which, barring an efficient relief valve 
such as a polar jet piercing the star, implies the outward 
transport of a large amount of energy. The energy liberated 
by the formation and initial growth of the black hole must 
be absorbed by the stellar envelope. While the black hole is 
still small — perhaps only a few hundred solar masses — the 
liberated energy becomes comparable to the binding energy 
of the envelope, which therefore inflates by a large factor. 
It probably does not disperse, since that would cut off the 
black hole growth and therefore the energy supply. 

This leads to a structure dubbed a 'quasistar' (Begel- 
man et al. 2006, 200 8), which resembles a r ed giant or 
Thorne-Zytkow object l|Thorne fc Zvtkowll 19771) except that 
the power source for inflating the envelope comes from black- 
hole accretion instead of nuclear shell burning around a com- 
pact core. As the envelope inflates, the core density declines, 
which reduces the black hole growth rate and the energy 
feedback rate until a stable equilibrium is reached when the 
energy output from black hole growth equals the Eddington 
limit for the total mass of the hole + envelope. The effective 
temperature of this bloated object is 

T QS ~ 5000e:i /5 m-^mI/ 6 20 K, (4) 

where the black hole and QS masses are expressed in units of 
10 4 Mq and 10 6 Mq, respectively. The parameter e — 0.1e_i 
represents an efficiency factor for the conversion of black hole 
growth into energy. In all of our models below, we adopt 
e_i = 1; if some of the energy were able to escape without 
percolating through the envelope, it would imply a smaller 
e and allow faster black hole growth during the QS phase. 
Since the black hole grows more rapidly than the QS 



envelope, Tqs decreases with time. But it cannot d r op be - 
low the temperature corresponding to the lHavashil (|l96lf ) 
track. For Population III abundances, the minimum temper- 
ature is about 4000 K (which we adopt as the fiducial value); 
for near-solar enrichment, this number drops to ~ 3000 K 
or slightly lower. Unlike gas pressure-supported red giants, 
which remain stable as they hover close to the minimum 
temperature, QSs lose dynamical equilibrium and disperse 
once they reach this limit. This is because QSs are supported 
against collapse by radiation pressure, which also governs 
internal energy transport — implying that there are not 
enough degrees of freedom available to allow the envelope 
to adjust quasistatically. From equation ([4]), we see that the 
black hole mass cannot exceed about 1 percent of the QS 
mass. For m > Se^T 44 / 9 , this leads to a final QS mass 
that substantially exceeds the mass of its supermassive star 
precursor, 

m.,8 ~ 1.4 e 4 _ /9 m 8/9 (5) 
and a black hole mass (at the time of dispersal) 

o -1/9 • 7/9 ra\ 

W!bh,4 ~ 3e_i m ' . (6) 

The lifetime of the QS phase is then 

iqs ~ lAe 4 _/?m- 1/g Myr. (7) 

For smaller infall rates, the final QS mass will be similar to 
the mass of the supermassive star (eq. [2]) and the black 
hole mass will be < 10 4 M . 



3 WHERE QUASISTARS FORM: MODELING 
THE HOST 

We investigate the formation and evolution of QSs and 
MBHs via cosmological realizations of the merger hierar- 
chy of dark matter halos from e arly times to the prese nt in 
a ACDM cosmology (WMAP5. Ikomatsu et aHl2009h . We 
track the merger history of 175 parent halos with present- 
day masses in th e range 10 11 < Mh < 1 15 M with a Monte 
Carlo algorithm (|Volonteri et al.ll2003l ) . The mass resolution 
of our algorithm reaches 10 5 Mq at z — 20, and the most 
massive halos are split into up to 600,000 progenitors. We 
also generate 12 Monte Carlo realizations of the merger hi- 
erarchy of a Mh = 2 x 10 13 Mq halo at z = 6, starting from 
z = 25. 

The formation of QSs requires large gas inflows, M > 
1 Mq yr~ , within the inner region (sub-parsec scales) of a 
galaxy. Following BVR2006 we associate these inflows with 
global dynamical instabilities driven by self-gravity, the so- 
called 'bars within bars' mechanism (jShlosman et al. i urn, 
1990). Self-gravitating gas clouds become bar-unstable when 
the level of rotational support surpasses a certain thresh- 
old. A bar can transport angular momentum outward on a 
dynamical timescale via gravitational and hydrodynamical 
torques, allowing the radius to shrink. Provided that the gas 
is able to cool, this shrinkage leads to even greater instabil- 
ity, on shorter timescales, and the process cascades. This 
mechanism works on a dynamical time (« Myr for high red- 
shift galaxies) and can operate over many decades of radius 
(from hundreds of pc to sub-parsec scales). 

Global bar-driven instabilities have now been observed 



4 Volonteri and Begelman 



in high-resolution numerical simulation s of gas-rich galax- 
ies dWise et alj|2008l: iLevine et al.ll2008l ; iRegan fc Haehneltl 
120091 : iMaver et alJl2009l ). While BVR2006 suggested that 
fragmentation and global star formation are suppressed only 
if the gas temperature remains close to the virial tempera- 
ture (e.g., due to lack of metals and molecular hydrogen), 
these simulations find strong inflows that occur before most 
of the gas fragments and forms stars also at solar metal- 
licities, possibly due to the continuous generat ion of super- 
sonic turbulence dBegelman fc Shlosmanll2009l ). In contrast 
to BVR2006, we therefore relax the assumption that QSs 
can be formed onl y out of zero-metallicity gas. Inspired by 
IMaver et all (|2009h . we instead assume that inflows are trig- 
gered by gas-rich major mergers. 

The gas content of galaxies is regulated by the competi- 
tion between stellar wind/supernova feedback that depletes 
galaxies of gas, and replenishment of the gas reservoir via 
mergers with gas-rich halos. In particular, we assume that 
Pop III stars form i n metal-free halos with T v i r > 2000 K 
ijYoshida et al.l l2006). where we model m etal enrichment by 
the 'h i gh feedback, best g u ess' m odel of IScannapieco et al.l 
( 2003). IScannap ieco et all (|2003l ) model metal enrichment 
via pair-instability supernovae winds, by following the ex- 
pansion of spherical outflows into the Hubble flow. They 
compute the comoving radius, at a given redshift, of an out- 
flow from a population of supernovae that exploded at an 
earlier time. Using a modification of the Press-Schechter 
technique (Scannapieco & Barkana 2002), they then com- 
pute the bivariate mass function of two halos of arbitrary 
mass and collapse redshift, initially separated by a given 
comoving distance. From this function they calculate the 
number density of supernovae host halos at a given comov- 
ing distance from a 'repicient' halo of a given mass Mh that 
form at a given redshift z. By integrating over this func- 
tion, one can calculate the probability that a halo of mass 
Mh forms from metal-free gas at a redshift z. When a halo 
forms in our merger tree we calculate the probability that it 
is metal-free (hence, it can form Pop III stars) and determine 
if this condition is satisfied using Monte Carlo techniques. If 
a Pop III star forms in a halo, we conservatively assume that 
all gas is expelled, leading to a nil gas fraction, / gas = 0. Gas 
is replenished only via mergers with gas-rich halos, that is, 
halos that have never experienced mass loss through Pop III 
outflows. When two halos merge, we sum their gas masses. 
We allow QS formation only in mergers where the final gas 
fraction exceeds a specified threshold, / gas > / t hr- We have 
varied /thr between 0.025 and 0.25 and found that the model 
results are largely insensitive to the choice of /thr- 

When selecting the sites of QS formation, we addition- 
ally require that the merger remnant has very low angular 
mome ntum, which ensures the optimal conditions for gas 
infall (Begelman et al. 20061'). We parameterize the angular 
momentum of a dark matter halo with mass Mh via its halo 
spin parameter, A sp i.n = JhE^ 2 /(GM^/ 2 ), where Jh is the 
total angular momentum and E^ is the binding energy. The 
exact angular momentum threshold below which inflows are 
triggered, A t hr, is a free parameter that we constrain against 
observations (see Section 5 below). 

iD'Onghia fc Navarrol |2007l ) investigate the spins of ha- 
los that are remnants of major mergers, and compare them 
to the 'global' population using cosmological N-body simu- 
lations. They find that halos that are still unrelaxed after a 



major merger tend to have higher-than-average spins. They 
quantify the spin parameter distributions for post-merger 
halos and relaxed halos. In both cases the distributions are 
log-normal, but with different parameters: Aspin = 0.028 and 
a\ — 0.58 for relaxed halos, while A sp m = 0.04 and ax = 0.65 
for unrelaxed systems. When a halo forms, we pick its spin 
parameter from the distribution for relaxed halos. After a 
major merger, we pick a new A sp m from the unrelaxed log- 
normal distribution, and compare it to A t hr- If the halo has 
a spin parameter below the threshold, we consider the halo 
a candidate for QS formation^- 

For halos that meet both criteria, / gas > / t hr and 
A sp in < Athr, we calculate the gas inflow rate as M = v't/G, 
where u c is the circular velocitjo If M > 1 Mq yr _1 , then we 
model supermassive star, QS and seed black hole formation 
as described in Section 2. We limit the mass of the form- 
ing QS to less than the mass of available gas, / gaa Mh. From 
stability considerations, we also assume that the supermas- 
sive star that precedes the QS undergoes core collapse if its 
mass exceeds 10 7 Mq even if it has not exhausted its core 
hydrogen. 

We further assume that if a seed MBH is already 
present, QS formation is suppressed. This criterion comes 
from the fact that any QS with a black hole mass exceeding 
about 1 per cent of the envelope mass will violate the min- 
imum temperature condition and be dispersed by radiation 
pressure (Begelman et al. 2008). Thus, it is very difficult 
to lay down a massive envelope of gas around an existing, 
naked MBH. 

Summarizing, our model assumptions are as follows: 

• inflow triggered by gas-rich major mergers (mass ratio 
between 1:10 and 1:1); 

• inflow if gas fraction / gas > / trir and spin parameter 
A ap in < Ath r ; 

• central mass accretion rate M = v% /G; 

• if M > 1 M /yr then QS and seed MBH form, with 
masses given by 

- Mqs = 1.4 x loV^m 8 / 9 M 

- M BH = 3 x 10 4 eli /9 m 7/9 M ; 

• QS formation suppressed if MBH already present in 
halo. 

We compare two models, one in which QSs are the only 
route to MBH seed formation, and one in which Pop III 
remnants offer a complementary channel. In both cases, we 

2 The spin parameter for the g aseous component, A gaB , c an be 
different from the A for the halo. iGottlober fc Yepesl (|2007| > find, 
however, that statistically A gaa /Ah a i ~ 1.4, and that the ten- 
dency is for this ratio to decrease toward unity as redshift in- 
creases. We expect, therefore, that statistically our approach is 
valid. 

3 A halo of mass Mh collapsing at redshift z has a circular velocity 



142 



I 1/3 



10 12 M 



1 1/6 



(l + 2 ) 1 / 2 kms- 1 (8) 



where A c is the over-density at virialization relative to the crit- 
ical density. For a WMAP5 cosmology we adopt here the fit- 
ting formula (Bryan & Norman 1998) A c = 18tt 2 + 82d - 39d? , 
where d = — 1 is evaluated at the collapse redshift, so that 

= n m (i + z) 3 /[n m (i + zf + + n k (i + z) 2 }. 
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assume that Pop III stars form in metal-free halos with 
T v h- > 2000 K0. In the QS-only scenario, we assume that Pop 
III stars form but leave no massive remnan t (that is, the star 
mass is below 260 Mm; [Fryer et all feoOll) . Therefore, Pop 
III stars in this model only provide stellar wind/supernova 
feedback that depletes galaxies of gas. In the Pop III rem- 
nants+QS scenario we assume a logarithmically flat ini- 
tial mass function, dN/d\ogM =const, be tween 10 Mp and 
600 M , where the upper limit comes from lOmukai fc Pallal 
|2003t ). and suppose that seed MBHs form when the progen- 
itor star is in the mass range 40 - 140 M or 260 - 600 M 
l|Frver et al.ll200ll ). The remnant mass is taken to be one- 
half the mass of the star. Once a seed MBH forms, we sub- 
sequently treat it in exactly the same way, regardless of its 
formation mechanism. 



4 HOW BLACK HOLE SEEDS GROW: 

MODELING THE COSMIC EVOLUTION 

We base our model for MBH mass growth on a set of sim- 
ple assumptions, suppor ted by both simula tions of quasar 
triggering and feedback (ISpringel et al.l f2005). and analysis 
of the relationship between MBH masses (Mbh) and the 
prop e rties of their hosts in today's galaxies jGebhardt et al.l 
|2000| ; iFerrarese fc Merrittll20od ; lFerrarese!l2002l ). The main 
features of the models have been discussed e lsewhere 
jVolonteri et al.ll2003l ; IVolonteri fc Nataraianll2009l . and ref- 
erences therein). We summarize in the following the rele- 
vant assumptions. MBHs in galaxies undergoing a major 
merger (i.e., having a mass ratio > 1 : 10) accrete mass 
and become active. Each MBH accretes an amount of mass, 
AM = 9 x lO 7 (cr/2OOkms~ 1 ) 4 ' 24 M , where a is the veloc- 
ity dispersion after the merger. This relation ship scales with 
the M bh — cr relation, as it is seen today l|Giiltekin et al.l 
120091 ): 



the normalization in AM was chosen to take into account 
the contribution of mergers, without exceeding the mass 
given by the Mbh — o relation. 

We link the correlation between the black hole mass 
and the central stellar velocity dispersion of the host with 
the empirical correlation between the central stellar velocity 
dispersion and the asymptot ic circular velocity ( v c ) of galax- 
ies (F errarese 2002; see also iPizzella et al.l 12005: Bac s et al.l 

Mm, 



a = 200 



( V S V 

V304 kms" 1 / 



, kms" 1 . (10) 
,304kms- 1 / v ' 

The latter is a measure of the total mass of the dark matter 
halo of the host galaxy. We calculate the circular velocity 
from the mass of the host halo and its redshift. 

4 We neglect the effect of a redshift-dependent Lyman-Werner 
background that progressively increases the minimum mass nec- 
essary for Pop III star formation (Trenti &; Stiavelli 2009). As 
shown below, the differences in the numbers of forming QSs and 
MBHs, for the cases with and without Pop III remnants, are neg- 
ligible. If the Lyman- Werner background reduces the number of 
forming Pop III stars, the result will fall in-between the two cases 
considered here. 



The rate at which mass is accreted scales with the 
Eddington rate for the MBH, and we set a fixed Edding- 
ton ratio of .fEdd = -3, based on observations of quasars 



l|Kollmeier et al.l 20061 ). This is also the strong upper limit 
for rad iatively efficient spherical accretion onto Pop III rem- 
nants i|Milosavlievic et all I2OO9T ). Accretion starts after a 
halo dynamical timescale and lasts until the MBH, of ini- 
tial mass Mi n , has accreted AM. The lifetime of an AGN 
therefore depends on how much mass it accretes during each 
episode: 



£agn = 



^Edd 



/Edd 1 



Tn(M fln /M in ), 



(11) 



where e is the radiative efficiency (e ~ 0.1), 
tedd = 0.45 Gyr and M fln = min[(M in + AM), 1.3 x 
10 8 (o-/200kms- 1 ) 424 M Q ]. 

We further assume that, when two galaxies hosting 
MBHs merge, the MBHs themselves merge within the 
merger timescale of the host halos, which is a plausi- 
ble assumption for MBH binari es formed after gas-rich 
galaxy mergers ijDotti et al.ll2007l . an d references therein) . 
We adopt the relations sugg ested by iTaffoni et all (|2003l ) 
for the merger timescale. 



5 MODEL RESULTS: QUASISTARS AND 
BLACK HOLES 

5.1 Constraining the quasistar and black hole 
formation efficiency 

As described in Section 3, our model depends mostly on one 
parameter, A t hr- We constrain A t h r from below by using the 
number density of 2 ~ 6 quasars. We obtain a (less certain) 
upper limit on Athr by estimating the maximum permitted 
contribution of supermassive stars to reionization, and later 
show (in Section 5.3) that this constraint is consistent with 
the Soltan (1982) constraint on integrated quasar light. 

The population of 2 ~ 6 quasars probes the high- 
luminosity end of the lumin osity function, ro ughly L 
10 ergs - . Recent surveys (jjiang et alj l2009t and refer- 
ences therein) yield a number density of these luminous 
sources, powered by billion-solar-mass MBHs (assuming 
sub-Eddington luminosities and negligible beaming), of ~1 
Gpc -1 . This comoving density roughly corresponds to the 
number density of halos with mass < 10 13 M in a WMAP5 
cosmology. With the scalings described in Section 4, we can 
express the mass of a halo hosting a 10 9 M MBH as follows: 

T m I - 84 To A 1 _1 / 2 

A/h = 5x10 [lond feis^J {1+z) Me ' 

(12) 

which is in basic agreement with the number density-based 
estimate. If Pop III remnants, with mass ~ 100 M , pro- 
vided the only route to MBH formation, equation (11) would 
imply that for an Eddington ratio of 0.3 and a radiative effi- 
ciency of 0.1, t agn = 2.68 Gyr is needed for a MBH to reach 
10 9 M . This is more than twice the age of the Universe at 
2 = 6. Even decreasing e to 0.05-0.06 (appropriate for ei- 
ther a Schwarzschild BH, or a BH fed by a counter-rotating 
accretion disc) does not bring the timescale within an ac- 
ceptable range (< 0.95 Gyr in a WMAP5 cosmology). For 
Pop III remnants to reach a billion solar masses by 2 = 6, 
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very specific conditions must be met: accretion occurs con- 
tinuously and at rates close to /ecm = 1, if not higher. If 
either of these conditions fails, Pop III remnants alone can- 
not reach the required MBH masses. Clearly, our approach 
here is very conservative. The rate at which MBHs in the 
early Universe accrete is unknown, and it is possible that the 
most massive high-redshift MBHs may be actually be fe d 
at the Eddington limit (or above: I Volonteri fc Reesll2005h . 
One open question is the role of feedback , if accretion is 
radiat i vely efficient. Pelupessv et al.l |2007l '): lAlvarez et al.l 
( 2009t); iMilosavlievic et al.1 |2009) suggest that the growth 
of the MBH is severely limited by thermal feedback to rates 
not higher than 30% of the Eddington rate. 

On the other hand, at very high accretion rates, the 
excess radiation might be effectively trapped, and a black 
hole accreting at super-critical values do es not necessar- 
ily radiate at super-Edd ington luminosity (Bcacrmanl ll979l ; 
iBegel man fe Meier! fl982l ). Also, when the accretion rate is 
high, the formation of collimated outflows is common (e.g. 
blazars). Such collimated jets may not cause feedback di- 
rectly on the host (which is pierced through), but deposit 
their kinetic energy at large distances , leaving the host un- 
scath ed (in a different context, see IVernaleo fc Reynolds! 
2006). This is likely if at large accretion rates photon 
trapping decreases the disk luminosity, while concurrently 
the presence of a jet helps to remove angular momen- 
tum, thus promoting efficient accretion (Volonteri, Ghis- 
ellini & Haardt, in preparation). This is supported by the 
strong cosmological evolution of blazars detected by the 
Burst Alert Telesc ope (BAT) onboard the Swift satellite 
ijAiello et al.l 120091 ) . The derived luminosity function indi- 
cates that beamed, jetted sources become dominant at the 
highest MBH masses (> 1O 9 M0) and redshifts (z > 5), when 
compared to the cos mic evolution of radio-quiet sources 
l|Ghisellini et al.ll2010l ). 

Figure [T] shows the mass growth of the central black 
hole in the main halo of merger trees corresponding to Mh = 
2 x 10 13 Mq halos at z = 6. We see that the requirement set 
by z = 6 quasars, powered by 10 9 Mq MBHs, requires A t hr 
to be larger than 0.01. 

The case A t hr = 0.01 is only marginally consistent with 
the constraints. This is because the probability of MBH for- 
mation in a given halo is less than 10 -2 , given the functional 
form of the spin parameter probability distribution. A halo 
needs ~ 100 major mergers to have a probability of MBH 
formation close to unity (provided, of course, that v c is large 
enough to provide the required M). In the 12 realizations 
of M h = 2 x 10 13 M halos at z = 6, we have not found a 
single case in which the main halo of the tree forms a MBH 
if Athr = 0.01. The MBH hosted in the central halo is in- 
stead acquired during a merger with another halo hosting a 
central MBH. This probability is lower if MBH formation is 
less efficient (lower A t h r ), leading to the main halo acquiring 
a MBH only at relatively late cosmic times (z ~ 12 — 14), 
thus postponing its growth. When Pop III stars leave behind 
massive remnants, the situation worsens. Since we assume 
that QS formation is suppressed in the presence of a MBH, 
even fewer halos can form massive seeds that can accelerate 
the growth of a MBH by boosting the merger rateQ z = 6 



We have here ignored the effect of 'gravitational recoils' on the 



quasars therefore provide a lower limit to QS and seed for- 
mation efficiency: A t hr must be larger than 0.01. 

We estimate the possible contribution of supermassive 
stars to hydrogen reionization as follows. Since the hydro- 
gen recombination timescale at z > 6 is shorter than the 
then-Hubble time, to derive an upper limit to the contribu- 
tion to reionization we can integrate the number of ionizing 
photons over time, that is, the sum all the ionizing photons 
emitted by all the forming supermassive stars. To estimate 
the number of ionizing photons emitted by a single super- 
massive star, we assume the photon output of a blackbody 
with temperature and radius given by T» and R, from Sec- 
tion 2, and multiply the output rate by the main sequence 
lifetime tMS (eq. [I]). Since the mass function of QSs — and 
the preceding supermassive stars — peaks at ~ 10 6 Mq (see 
Section 5.2 below), we calculate the photon output for rh = 1 
and use T c> g = 1.4 as before. This gives an ionizing photon 
output of 1.6 x 10 67 ionizing photons per supermassive star. 

The total number density of ionizing photons emitted 
by the supermassive star population can then be approxi- 
mated as: 



«io 



(z) « 1.6 x 10 67 J2 n SMS,i ph Mpc" 



(13) 



i=0,j{z) 



where the sum is over the number density of supermassive 
stars (nsMs) that form during each timestep between our 
starting redshift (z max = 20, timestep 0), and the timestep 
j(z) corresponding to the redshift of interest, z. 

We compare this number density of photons to the num- 
ber density of baryons: 



/bar ^7 



Pcrit 

(m) 



7.4 x 10 66 Mpc" 3 , 



(14) 



where /bar — 0.17, S7„ 
and (m) ~ 1.7 x 10~ 24 g. Therefore, 



0.3, pcrit = 9.2 x 10" 30 g cm" 3 , 



nion(z) 
^bar 



1.6 
0.74 



nsMS,i- 



(15) 

This calculation in shown in Figure [2] If we consider a reion- 
izatio n criterion 1 < ni on /nb ar < 3 at z = 6 (|Gnedinl 
2008), we see that A t hr > 0.02 tends to overproduce ion- 
izations, provided that the ionizing radiation escapes from 
the halo core. As noted in Section 2, this is far from cer- 
tain because the gas falling onto the supermassive star is 
optically thick. Our estimate of the ionization rate should 
therefore be regarded as an upper limit, corresponding to 
the case where the photospheric radiation of the supermas- 
sive star is unobscured (e.g., as the result of disc- like inflow) 
or at least not degraded. The ionization rate from super- 
massive stars therefore provides a weak constraint. However, 
we note that an additional limit on QS and BH formation, 



MBH population. For MBHs merging in a gas-rich environment , 
as expected for these high-redshift systems, Volo nteri et al.l [120101 ) 
have shown that the probability of ejection is negligible. This is 
because the accreting gas exerts gravitomagnetic torques that suf- 
fice to align the spins of both MBHs with the angular momentum 
of the large-scale gas flow, leading to aligned spin-orbit configura- 
tions that yield recoil velocities < 100 km s — 1 . Even for the 'worst 
case scenario' of isotropically distributed spins, the probability of 
ejection for MBH binaries in the massive haloes which host very 
high redshift quasars effectively drops to zero at z ~ 13 — 19. 
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discussed in Section 5.3 (Soltan's argument), reinforces the 
limits discussed in this section. In the following we assume 
0.01 < A t hr < 0.02, and show results for Athr = 0.01 (low 
efficiency) and A thr = 0.02 (high efficiency). 

The quasar population that results from MBHs formed 
in QSs does not contribute significantly to hydrogen reion- 
ization (see Fig. [3]) for the low efficiency case, while the 
high efficiency case produces many more MBHs that pro- 
duce significant ionizing radiation. High efficiency MBH for- 
mation (Athr = 0.02) is also constrained to be an upper 
limit by the contribution of quasars to helium reionization. 
Taking into consideration ionization energy thresholds for 
Helfl and H (54.4 vs. 13.6 eV) and number densities of the 
two species (12 times more H than He), there are roughly 
3 times more He ionizing photons per He atom than pho- 
tons above 13.6 eV per H atom for our quasar spectrum. 
As He recombines ~6 times faster than H, the number of 
ionizing photons per He atom emitted in one recombina- 
tion time is only roughly half of the corresponding value 
for H, and there will be just a small delay between H and 
He reionizations. In the low-efficiency model He reioniza- 
tion occurs at z ~ 3, while in the high-efficiency case He is 
reionized at z ~ 3.5 — 5 (depending on the dumpiness of 
the intergalactic medium and the escape fraction of ionizing 
photons). Although observational constraints on He reion- 
ization are still wea k, they suggest that Hell r eionization 
occurs at z ~ 3 — 4 (Fa ucher-Giguere et alJl200St and refer- 
ences therein). Hence, while the low-efficiency model is con- 
sistent with He reionization, the high-efficiency case clearly 
represents a strong upper limit to MBH seed production, 
as, if MBHs accrete as efficiently as assumed here, He is 
reionized at much higher redshift. We note that hard pho- 
tons from supermassive stars do not contribute much to He 
reionization (the number of ionizing photons for He atom is 
roughly one-third of the number of ionizing photons for H 
atom shown in Fig. • 

5.2 Mass functions and formation times of 
quasistars 

Figure [3] shows the mass function of QSs, integrated over 
all redshifts of formation. The mass function peaks at ~ 
10 6 Mq , corresponding to the minimum gas infall rate M = 
Vc/G = IMq yr _1 . Since the infall rate is related to the 
circular velocity of the host halo, it is also linked to the mass 
of the halo. The halo mass function in a ACDM cosmology 
is a steep function of the halo mass; the QS mass function 
therefore peaks at the lowest allowed mass (~ 10 6 M ). 

A QS ends its life with a photospheric temperature 
of about 4000 K for Population III abundances (which we 
adopt as the fiducial value). We model the emission as a 
blackbody at Tqs = 4000 K for the entire lifetime of the 
QS (equation 7), and calculate number counts in the 2-10 
|im band of the James Webb Space Telescope (JWST), as- 
suming a sensitivity of 10 njy at 2/im. These are shown in 
Figure[4]for our reference models. In the high-efficiency case 
(Athr = 0.02), JWST could detect up to a few QSs per field 
at z ~ 5 — 10. In each panel of Figure 01 the lower curves 

6 The ionization threshold for He I, 24.6 eV, is low enough that 
He I and H reionization are expected to occur almost concurrently. 




10 15 10 15 

redshift redshift 



Figure 1. Growth of the central black hole in the main halo of 
merger trees corresponding to Mh = 2 X 10 13 Mq halos at z = 6 
(we show here only 5 realizations out of 12 for clarity). Lower 
panels: relatively low-efficiency QS/supermassive star formation 
case (A t hr = 0.01). Upper panels: efficient QS/supermassive star 
formation case (A t hr = 0.02). Left panels: Pop III stars leave no 
remnant that can impede QS and MBH seed formation. Right 
panels: Pop III stars leave remnant MBHs. The big circle is cen- 
tered at z=6 and M B h = 10 9 M Q . 



show the expected number counts for fields at the depth 
of the Hubble Ultra Deep Field, observed with WFC3 on 
board of the Hubble Space Telecope in the F160W band, to 
an AB magnitude limit of 28.8 (Bouwens et al. 2010). We 
conclude that it is possible that QSs can already be seen in 
deep WFC3 images, although it is not very likely. 

Although the peak of QS formation occurs at z > 5, 
some QSs form at redshift as low as z ~ 2. These late 
bloomers form in halos with relatively low cosmic bias. We 
recall that the bias of a halo is the number of standard de- 
viations that the critical collapse overdensity represents on 
the halo mass scale, Mh- v Q = 8 c (z)/a(Mi 1 , z). At any red- 
shift we can identify the characteristic collapsing mass (i.e., 
v c = 1), and its multiples. The higher v c , the more massive 
and rarer the halo, and the higher its bias and clustering 
strength. Figure [5] exemplifies our results for two cases: a 
Milky Way-sized halo (M h = 2 x 10 12 M Q at z = 0, lower 
bias), and a cluster-size halo (M h = 1O 15 M at z = 0, 
higher bias). We therefore expect that the latest-forming 
QSs should be found in the field, rather than in a high- 
density environment. As one goes to lower and lower galaxy 
masses, the number of halos where QSs can form decreases, 
as the requirement M oc v% > 1 Mq jyr selects only halos 
with mass larger than ~ 10 9 Mq as possible QS/MBH hosts. 
Given that the probability of QS/MBH formation in a given 
halo is 1-10% (due to the log- normal distribution of spin 
parameters), and that major mergers between ~ 10 9 Mq 
halos are very rare events in the merger history of galax- 
ies with masses below ~ 10 11 Mq, we expect dwarf galaxies 
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Figure 2. Number of reionizing photons per baryon for two 
QS formation efficiencies, A t i lr = 0.02 and A t j lr = 0.01. Cases 
where Pop III remnants offer an additional channel for MBH 
formation yield almost identical results. Solid curves: supermas- 
sive star contribution to the reionizing budget. Dashed curves: 
contribution from quasars associated with QS remnants where 
we assume a classic m ulticolor disc spectrum up to fcT ma x ~ 
1 keV (Mbh/ Mq) -1 / 4 iShakura fc Sunvaevl Il973h . and a non- 
thermal power-law component L„ oc v~ a , with a fs 1 at higher 
energies. With this choice, a MBH with mass ~ 10 4 Mq ra- 
diates 20% of its bolometric power as hydrogen-ionizing pho- 
ton s, with (E~f) = 140 eV ( for details on the implementation 
sec I Volonteri fc Gnedinl [200j) . Dotted curve: contribution from 
normal stars (from Gncdin 2008). 

to be rarely seeded with MBHs (results for the frequency 
of MBHs in dwarf galaxies are very similar to those found 
in I Van Wassenhove et ai1l2010t ). In order to check whether 
the infrequency of QS/MBH formation at low halo masses 
is compensated by the larger abundance of small galaxies, 
we performed 150 merger trees each for halos with mass 
5 x 10 9 M , 10 10 M and 5 x 10 10 Mq. We find that forma- 
tion of QS/MBHs becomes statistically negligible for halos 
with mass below 10 10 Mq , and that QSs in low- mass galaxies 
do not dominate the overall population. 

Figure [3] also shows the mass function of seed black 
holes, integrated over all redshifts of formation. The mass 
function peaks at ~ few 1O 4 M , but rare supermassive 
seeds, with masses ~ 1O 6 M , are possible. The more ef- 
ficient QS and BH formation is, the steeper the mass func- 
tion. This is because more seeds form at the highest redshifts 
when the hosts have relatively low mass and low circular ve- 
locity 

5.3 Integrated properties of the black hole 
population 

Figure [6] shows the number density of MBHs: depending on 
the efficiency of QS and MBH formation (A t hr), the number 
density can be dominated by either Pop III or QS MBHs. In 



Figure 3. Left: mass function of QSs, integrated over all redshifts 
of formation. Right: mass function of seed black holes. Lower set 
of histograms: relatively low-efficiency QS formation case (A t j lr = 
0.01). Upper set of histograms: efficient QS formation case (Ay lr = 
0.02). Dashed curves: Pop III stars leave no remnant that can 
impede QS and MBH seed formation. Solid curves: Pop III stars 
leave remnant MBHs. 

the low-efficiency case (Athr = 0.01) the number density of 
MBHs is dominated by the Pop III channel, if available; the 
opposite is true for the high-efficiency case. Note how the 
light Pop III remnant seeds are much more often actively 
accreting (compare left and right panels). This is an effect 
of imposing self-regulation of the MBH mass, i.e., that in 
a given accretion epis ode a MBH cannot grow to exceed 
the Mbh — a relation (|Volonteri fc Nataraianl [20091 ). Since 
QS MBHs are born above the expected Mbh — a relation 
for their host, accretion is effectively damped. Light Pop III 
remnants, instead, are undermassive and accrete steadily. 

The number density of MBHs at z — can be compared 
to the number density of galaxies, using as a constraint, 
for instance, the galaxy luminosity function as measured in 
the Sloan Digital Sky Survey (e.g., Benson et al. 2007). At 
stellar masses ~ 10 8 M the galaxy number density is ~ 
0.3— lMpc~ 3 , which is well above our MBH density at z = 
(see also Figure 11 in Carlberg et al. 2009). 

In contrast to the number density of MBHs, the mass 
density is dominated by QS MBHs (Figure [7J). We show 
both the mass density locked into MBHs and the cumula- 
tive mass d ensity accret ed onto MBHs (related to the Soltan 
constraint: ISoltanlll982b . Soltan's argument offers an addi- 
tional constraint on our models. We compare the cumulative 
ma ss density accreted o nto MBHs with that derived from 
the lHopk ins et all (|2007l ) bolometric luminosity function of 
quasars: 

PQSo(< Z) = 5-^ / / , dl 

« 2 J J h (i + z)y/n m (i + zf + n A 

(16) 

where the mass accretion rate, M acc = Le 1 c 2 , is con- 
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Figure 4. Number counts of QSs for a JWST sensitivity of 10 
njy at 2/xm (upper curves) and for the HUDF WFC3/IR in the 
F160W band, to an AB magnitude limit of 28.8 (Bouwens et al. 
2010: lower curves). We assume Tqs = 4000 K. Lower panels: 
relatively low-efficiency QS formation case (A t hr = 0.01). Upper 
panels: efficient QS formation case {X^r — 0.02). Left panels: Pop 
III stars leave no remnant that can impede QS and MBH seed 
formation. Right panels: Pop III stars leave remnant MBHs. The 
dotted horizontal line shows the limit of 1 source per JWST field 
(2.2 X 2.2 arcmin 2 ). 

_ i i i | i i i | i i i | i i i | i i i | i i i | i i i | i _ 
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Figure 5. Formation times of QSs for two differently biased sys- 
tems: a MW-size halo and a cluster-size halo. Dot-dashed curves: 
Pop III stars leave no remnant that can impede QS and MBH 
seed formation. Solid curves: Pop III stars leave remnant MBHs. 
Here TVtot is normalized to a comoving cubic Mpc. 



verted into MBH mass growth, Mbh = (1 — e)Le _1 c -2 , with 
e the energy conversion coefficient. The boundaries assume 
a minimum radiative efficiency e = 0.06 and a maximum 
radiative efficiency e = 0.3 (Schwarzschild to rapidly rotat- 
ing Kerr BH). All our models fit comfortably within the 
allowed range. We see from Figure [7]that Soltan's argument 
provides an additional constraint to the upper limit on MBH 
formation via QS seeds. If the QS formation efficiency were 
much larger than A t hr = 0.02, the cumulative accreted mass 
density at z — 4 would be overestimated. 

Note that the value of the mass density locked into 
MBHs at z = 0, p BH .q = (3.2 - 5.4) x 10 5 M Q Mpc -3 (e.g., 
lMerloni~fe Heinz 200a an d references therein) is another, in- 
dependent constraint that can be evaluated from Figure [7] 
The low-efficiency case is only marginally consistent with 
Pbh,o- For comparison we show here a case of MBH evolu- 
tion that considers the Pop III remnant channel only. This 
is also hardly consistent with pbh,o and Soltan's argument. 
If BH seeds are rare or very light, a constant /Edd = 0.3 is 
insufficient to explain the growth of the BH population and 
its current demography. 

The mass function of MBHs is show n in Figure [8] We 
also s how the mass function predicted bv lMerloni fc Heinz! 
(2008 ), derived using a completely different technique. 
iMerloni fc Heinj l|2008f) solve the continuity equation for 
the MBH mass function using the z = one as a bound- 
ary condition, and the luminosity function as tracer of the 
MBH growth rate. There is no a priori reason to expect that 
enforcing the A/bh~ a relation and assuming a fixed Edding- 
ton rate would match the z > m ass function of MBHs pre- 
dicted bv lMerloni fc Heinj (|2008T ). The very good agreement 
at masses > 10 7 M© suggests that, albeit very simple, our 
accretion prescription is well motivated. The disagreement 
at low masses could be due to incompleteness in the Merloni- 
Heinz inventory of low-mass MBHs, since their model works 
backwards from z — 0, and it does not assume an initial 
seed mass below which a MBH mass cannot decrease. If the 
only channel of MBH formation is via QS seeds, then the 
mass function cuts off below log^A/sH/M©) = 4.5, where 
the seed mass function drops as well. If Pop III remnants 
offer an alternative route to MBH formation, then the mass 
function is double-peaked, with each peak corresponding to 
a different mechanism for the formation of seeds. The peak 
at log(A/aH/M0) ~ 4.5 becomes more pronounced with a 
higher efficiency of QS (and resulting seed BH) formation. 
Additional factors, such as mass and redshift dependent ac- 
cretion rates (we used a simple fixed-Eddington rate model) 
and obscuration, will modulate the exact shape of the mass 
function of observable MBHs. 



6 DISCUSSION AND CONCLUSIONS 

We have used a merger-tree approach to estimate the rate at 
which supermassive (> 10 6 M ) stars might have formed as 
a function of redshift, and the statistical properties of the re- 
sulting quasistars and seed black holes. Key to the formation 
of supermassive stars is a large gas infall rate (M > 1 Mq 
yr _1 ), driven by global gravitational instability in the poten- 
tial of a recently merged halo. We use the analyses presented 
by Begelman et al. (2008) and Begelman (2010) to estimate 
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Figure 6. Number density of MBHs vs. redshift, integrated over 
all MBH masses. Lower panel: relatively low-efficiency QS forma- 
tion case (A thr = 0.01). Upper panel: efficient QS formation case 
(Athr = 0.02). Left panels: Pop III stars leave no remnant that 
can impede QS and MBH seed formation. Right panels: Pop III 
stars leave remnant MBHs. Thick curves (the lowest curve in each 
set): active MBHs only. Note that we do not include information 
regarding lifetime or a luminosity threshold. Depending on the 
efficiency of QS and MBH formation (Athr)i the number densities 
can be dominated by Pop III or QS MBHs (e.g., QS MBHs always 
dominate in the A t h r = 0.02 case, while in the A t h r = 0.01 case 
Pop III remnants can dominate the MBH population in number.) 



the properties of the resulting supermassive star, QS and 
seed black hole as a function of this rate. 

We relate the properties of the seed black holes to the 
observable distributions of active and quiescent black holes 
at different redshifts, by applying a simple set of rules for 
their subsequent growth by accretion and mergers. We com- 
pare our direct collapse seed model to a pure Pop III seed 
model, and also study the case in which both channels for 
MBH formation operate simultaneously. 

Our results can be summarized as follows. 

• A limited range of supermassive star/QS/MBH forma- 
tion efficiencies exists that allows one to reproduce observa- 
tional constraints (the density of z = 6 quasars, the cumu- 
lative mass density accreted onto MBHs, the current mass 
density of MBHs, reionization) . These constraints translate 
into 0.01 < A thr < 0.02. 

• The mass function of QSs peaks at Mqs ~ 10 6 M© , and 
decreases almost as a power-law with slope ~ —(1.5 — 2). 
The more efficient QS and MBH formation is, the steeper 
the mass function. This is because more seeds form at the 
highest redshifts when the hosts have relatively low mass 
and circular velocity. 

• Modeling QS emission as a blackbody at Tqs = 4000 
K, we calculate the number counts for the JWST 2—10 fim 
band. Assuming a sensitivity of 10 njy at 2^tm, we find that 
JWST could detect up to several QSs per field at z ~ 5— 10. 



Figure 7. Mass density of MBHs vs. redshift. Panels and lines 
as in Figure [6] The thick line in each panel shows the integral 
over time of the accreted mass (cumulative accreted mass den- 
sity). The dashed curve in the upper-right panel shows the total 
mass density in MBHs when we consider Pop III remnants only. 
Hatched area: Soltan's constraint from the luminosity function 
of quasars using the Hopkins et al. (2007) bolometric luminosity 
function. The boundaries assume a minimum radiative efficiency 
e = 0.06 and a maximum e = 0.3 (Schwarzschild to rapidly rotat- 
ing Kerr MBH). 

• The redshift of formation of QSs increases with the cos- 
mic bias of their hosts. We therefore expect that the latest- 
forming QSs should be found in the field, rather than in a 
high-density environment. 

• The number density of MBHs can be dominated by ei- 
ther the descendants of Pop III remnants or QS MBHs. In 
the low-efficiency case (A t hr = 0.01) the number density of 
MBHs is dominated by the Pop III channel, if available; the 
opposite is true for the high-efficiency case (Athr = 0.02). In 
contrast to the number density of MBHs, the mass density 
is always dominated by QS MBHs. 

• If the only channel of MBH formation is via QS 
seeds, then the mass function of MBHs cuts off below 
^(Mbh/Mq) = 4.5, where the seed mass function drops 
as well. If Pop III remnants offer an alternative route to 
MBH formation, then the mass function is double-peaked, 
with each peak tracing a different seed formation mecha- 
nism. The peak at ^(Mbh/Mq) ~ 4.5 becomes more pro- 
nounced, the more efficient QS formation is. 

This model differs from previous proposed mechanisms 
(e.g. Lodato & Natarajan 2006) in various respects. First, 
gas accumulation in the central regions of protogalaxies is 
described by global, rather than local, instabilities (e.g., via 
the Toomre stability criterion). This implies that the spin 
parameter threshold for collapse cannot be derived by de- 
manding that the Toomre parameter, Q = (where £ is 
the surface mass density, c s is the sound speed, k = %/2Vh/R 
is the epicyclic frequency, and V is the circular velocity of 
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Figure 8. Mass function of MBHs in three redshift bins: 0.95 < 
z < 1.025; 1.9 < z < 2.1 and 2.8 < z < 3.2. The bin width 
was chosen to span similar cosmic times (~ 0.32 Gyr). In each 
panel we show two curves, the upper curve is the mass function 
of all black holes. The lower curve is the mass function of active 
black holes. We also show (big dots) the mass function predicted 
by Merloni & Heinz (2008). The very good agreement at masses 
> 10 7 Mq suggests that our MBH growth prescription is well 
motivated. Thick lines in the bottom panels: IMF of seed MBHs 
left behind by QSs. Left: A thr = 0.01 + Pop III MBHs. Middle: 
A thr = 0.02 + Pop III MBHs. If we switch off the Pop III channel, 
the mass function drops to zero at log(MeH/M0) = 4.5, where 
the seed mass function drops as well. Right: Pop III MBHs only. 

the disc), approaches a critical value, Q c , of order unity. 
We instead assume a fixed spin parameter threshold (cf. 
BVR2006). Second, we have here relaxed the assumption 
that large gas infall rates can occur only at zero metallicity 
to avoid fragmentation. The predicted mass functions have 
distinctly different shapes (compare Figure 3 with Figure 2 
in Lodato & Natarajan 2007), and MBH seeds tend to form 
later in the model presented here (z ~ 5 — 10 instead of 
z ~ 16 — 18, when zero metallicity is required). Differences 
between models can possibly be looked for at the highest 
redshifts, as at later times the growth of MBHs by accretion 
and mergers likely washes out dissimilarities. 

The most direct prediction of this work is shown in Fig- 
ure^ where we estimate the detectability of QSs in a JWST 
field of view. At low metallicities, QSs at redshifts of a few 
may resemble featureless blackbodies, with colors reminis- 
cent of brown dwarfs. Direct redshift measurements may not 
be feasible, in which case quasistars would have to be iden- 
tified via their massive hosts. Their low numbers — a con- 
sequence of their short lifetimes — will make it even more 
challenging to find them with a telescope having a relatively 
small field of view. It is possible that some QSs could have 
formed in metal-enriched regions, and even at relatively low 
redshifts (z > 2). These might be relatively easy to detect, 
but extremely rare. Other, secondary characteristics of QSs 
might conceivably aid in their detection, e.g., if accretion 



onto the black hole deep in the core leads to the formation 
of a jet that pierces the stellar surface. 

The shape of the MBH mass function below 
log(A/BH/Mo) = 4.5 is an additional diagnostic. Direct dy- 
namical measurement of MBH masses at the low-mass end 
is extremely difficult, due to the necessity of resolving the 
dynamics of stars within the sphere of influence of the puta- 
tive MBH, i.e., the region where the Keplerian potential of 
the MBH dominates over the overall galaxy potential. This 
region typically subtends much less than an arcsecond in 
nearby galaxies, at or below the resolution limit of existing 
8-10 m class telescopes. Future 20-30 m telescopes are likely 
to increase the sample of low-mass MBHs. 

Gravitational waves produced during the inspirals of 
compact objects into MBHs — extreme-mass-ratio inspirals 
(EMRIs), in particular — are expected to provide accurate 
constraints on the mass function of black holes at low red- 
shift. The proposed space-based gravitational wave detector, 
the Laser Interferometer Space Antenna (LISA), can probe 
the m ass function in the lO 4 A/0-lO 7 Af0 range. iGair et aD 
(2009) find that with as few as 10 events, LISA should con- 
strain the slope of the mass function below ~ 10 6 Mq to a 
precision ~ 0.3, which is the current level of observational 
uncertain ty in the low-mass slope of the black hole mass 
function l|Greene fc Holl2007r ). The combination of electro- 
magnetic and gravitational wave observations in the coming 
years will improve the currently limited constraints on what 
route, or routes, lead to MBH s eed formation. 

Numerical simulations (e.g.. lMaver et aDl2009l ) can test 
how strongly the 'bars within bars' inflow mechanism is tied 
to the angular momentum of the merging galaxies, and ei- 
ther validate our hypothesis that a single parameter A t i lr 
can describe the efficiency of the cascade process, or indi- 
cate that additional parameters influence the formation of 
quasistars and their nested black holes. Even the existence 
of a second parameter may not have a strong effect on the 
results: for instance, we found that our results are not very 
sensitive to varying the threshold gas fraction, as long as 
it is above / gas ~ 0.025. It will also be necessary to study 
the global dynamical behavior of self-gravitating inflows in 
the inner regions, where the total gravitational potential is 
expected to approach a Keplerian shape and the dominant 
means of angular momentum transport may change from 
one involving large-scale bars (m = 2) to a flow resembling 
an eccentric disk (m = 1: e.g., Hopkins and Quataert 2009). 

It is even more critical to check how much of the in- 
falling gas actually reaches the small radial scales neces- 
sary to build the supermassive star. We have suggested that 
fragmentation and star formation may be less important 
than previously thought in quenching such inflows. Testing 
these hypotheses will require numerical experiments with 
high spatial resolution and a large dynamic range. 
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